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ABSTRACT 
A list of seismic events for the period 20 February to 19 March 1972 (The 
International Seismic Month - ISM) has been compiled using short-period data from a 
large number of seismic stations and arrays. Hypocenters and other parameters have 
been obtained for every event and are presented along with comparisons with other 
event lists for the same time period. Also included are the results of experiments 
designed to calibrate certain parameters which have been selected to indicate the 


accuracy Of cstimated depths and epicenters. The ISM list contains 996 verified events. 
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INTERNATIONAL SEISMIC MONTH EVENT LIST 


by R. T. Lacoss, R. E. Needham and B. R. Julian 


Introduction 


In January 1972 a group of seismologists from several countries met in Cam- 
bridge, Massachusetts,to discuss problems related to the seismic detection and identi- 
fication of underground nuclear tests. Considerable interest was expressed at that 
time in a proposal that a cooperative attempt be made to better assess the capabilities 
of deployed seismic instrumentation for monitoring worldwide seismic activity. Shortly 
thereafter several groups of seismologists agreed to informally collaborate in such an 
undertaking. The period from 20 February to 19 March 1972 (The International Scismic 
Month - ISM) was selected for the experiment and the Lincoln Laboratory agreed to act 
as a data center. 

The cooperating groups, more background, and some of the details of the Lincoln 
(1) 


role can be found in an earlier report Other reports of Lincoln activity and progress 


(2), (3), (4). 


can be found in three Semiannual Technical Summaries In the present 
report we give the complete and final list of events which occurred during the ISM 

and which we have been able to detect and locate using available seismic data. Also 
included for each event are the parameters assigned to the same events by the USGS 
PDE (U. S. Geologic Survey Preliminary Determination of Epicenters), the LASA (Large 
Aperture Seismic Array in Montana) bulletin, the NORSAR (Norwegian Seismic Array) 
bulletin, and an bulletin generated by a selected set of only 32 stations. Comparison 


between the ISM list and these others are given. No long-period data is presented or 


analyzed in this report. 


The Data 

Appendix I contains the list of 996 ISM events. The data for each event are listed 
in columns. For each event there are 60 parameters given. These parameters include 
the actual ISM parameters, USGS PDE parameters, LASA and NORSAR bulletin para- 
meters, and parameters obtained using a selected set of 32 stations. Parameter names 
are listed in the leftmost column of each page and the definition of cach of the parameters 
is presented at the start of the Appendix. Figure 1 is a map of the 996 epicenters in 
the ISM event list. 

Appendix II is a chronological listing of all LASA and NORSAR bulletin entries 
which could not be verified or associated with any of the 996 ISM events. These events 
may be real but we have not been able to verify them with data from other stations. 
Table | shows in summary form the number of events from various event lists which 


correspond to ISM events and which are unverified by the ISM list. 


TABLE I 

BIN TRIES IDENTIFIED UNVERILE TED 
SOURCE WITH ISM EVENTS ENTRIES 
USGS PDE 3953 Z 
LASA, BULLETIN ay 96 
NORSAR BULLETIN 361 oi) 
32 STATION ISM 818 () 
ISM LIST 996 0 


One event located by USGS using no station more than 3° from epicenter was excluded 
from the ISM as purely local. The other, in South Africa, required usc of Ps for location 


and our program did not have this option. 


Almost 200 stations were used in generating the ISM event list. When the list 
was completed the contributions of all stations were reviewed and a set of only 32 
was selected to determine the utility of such a smaller network. Station selection was 
ad hoc to give good global coverage while using stations which contributed data to large 
numbers of ISM events. The list of 32 stations and some indication of their contribution 
to the full ISM are given in Table Il. The rightmost column lists the number of ISM 
events with which the arrivals from that station could be associated. The next column 
tells how many arrivals were associated with the ISM events. The first numerical 
column gives the number of arrivals reported by the station. In the case of NORSAR 
and LASA this last is very large because reports of signal arrivals were obtained from 
internal array detection logs as well as from the array bulletins. In this report any 
reference to a 32 station list is to the list of events which could be generated using only 
data available from these 32 stations. 

Of all of the short-period stations used in the ISM only LAO, NAO, HES, KBL, 
and CHG and all Canadian stations were specifically read by an analyst or in some 
other way treated especially for the ISM. For all other stations only picks routinely 
submitted by the station to USGS for PDE were used. In general, as can be seen to some 
extent in Table II, the stations given special treatment tended to contribute to more 
events and to contribute more arrivals to the data base. We would have wished that 
all stations could be carefully picked for small events and that recording quality could 
be improved and filtering could be done. Unfortunately this was not possible and our list 


represents Only the best we could do with available data and resources. 


TABLE II 


DISTINCT ASSOCIATED ASSOCIATED 
STATION LOCATION INPUT PICKS PICKS EVENTS 
LAO Montana 3881 1358 734 
YKA N.W. Terr. Canada 634 509 500 
UBO Utah 1082 480 466 
NAO Norway 2297 582 426 
HFS Sweden 502 300 341 
MBC N.W. Terr. Canada 442 332 331 
KBL Afghanistan 983 298 292 
ASP N. Terr. Australia 555 281 281 
MAT Honshu, Japan 544 230 ye 
COL Alaska 318 231 228 
CHG Thailand 892 256 222 
PNS Bolivia 390 173 173 
ClA Queensland, Australia 304 172 171 
BLC N.W. Terr. Canada 196 167 166 
NUR Finland 215 166 162 
LUG Arizona LI? 152 150 
FBC N.W. Terr. Canada 166 147 147 
GBA India 268 144 144 
SPA South Pole 192 141 141 
PMG New Guinea 236 LZ 125 
KIC Ivory Coast Zot 114 112 
CLE E. Germany 237 109 107 
SSF France 123 103 100 
BDF Brazil 272 101 101 
BNG Cen. African Rep. 209 99 99 
QUE Pakistan 14] 99 99 
SHI Iran 185 79 78 
BAG Philippine Is. 84 69 68 
BUL S. Africa 118 68 66 
AFI Samoa 166 56 56 
EZN Turkey 144 2° 59 
SH. India 114 52 52 


Event Grades 

The origin times, hypocenters, and My, values for events in the ISM are not all 
equally reliable estimates of those parameters for different events. An attempt has 
been made to assign grades to events according to the probable reliability of the listed 
ISM event parameters. Eight grades have been defined and are included in the data of 
Appendix I. Grades A, B, C, D are for events which could be located without arbitrarily 
restraining depth. Whenever possible each ISM event has been located by an unrestrained 
nonlinear weighted least squares fit of reported arrival times and horizontal phase 
velocities at arrays. All major phases, including depth phases such as pP and pPKP, 
were treated in the same way by the programs. Any event for which convergence could 
be obtained with depth greater than zero and less than seven hundred kilometers was 
assigned a grade of A, B, C, or D. Events which could be located only by arbitrarily 
restraining depth to zero, seven hundred, or some value in between was assigned a 
grade of Al, BI, CI, or DI. The I indicates that the depth is indeterminate in the sensc 
just described. 


The event grades are based upon six different parameters. These are: 


NAVE MB: The number of individual station mM), values averaged 
to obtain the ISM my, value. 
DE GE : Degrees of freedom in the location and origin time 


estimates. This is equal to the number of scalar 
observations used in the location minus the number 
of parameters estimated (three for Al, BI, CI, Dl 
events and four for A, B, C, D events). 


RATIO: A measure of the goodness of fit of the observations to 


the estimated origin time and location parameters. 
Residuals were weighted inversely by their assigned 
standard deviations, squared, summed, and divided by 
DEGF to obtain RATIO. 

MAXAX2: Largest principal semiaxis of our estimated epicenter 
error ellipse. This parameter has kilometers as units and 
can be calculated only for A, B, C, or Devents. Results 
of some experiments designed to calibrate this parameter, 
and DEPTHQ which is discussed below, are presented in 
the final section of this report but more research is re- 
quired to establish its true meaning and utility. Additional 
discussion of the error estimation problem will be found 
in Reference 3. 

DEPT HO: Half length of our estimated depth error bar. It is cal- 
culated only for A, B, C, D events and has units of kilo- 
meters. It indicates the relative reliability of depths but 
like MAXAX2 requires further research for absolute cal- 
ibration and verification of its ultimate utility. 

DEPTH: The ISM assigned depth for an event in kilometers. 

The exact definition of grades in terms of the above parameters are given in Table III. 
For each grade all conditions must be met. The number of events in each grade is 


given in parenthesis in the table. 


TABLE J 


Grading Rules 


Depth Free for Convergence Depth Fixed for Convergence 
Grade = A (306) Grade = Al (71) 
MAXAX2 < 100 RATIO s 1 
RATIO =2 DEGF 25 
DEGF 24 NAVEMB 23 
DEPTHQ/DEPTH=0.5 or DEPTHQ=25 
NAVEMB 23 
Grade = B (34) Grade = BI (16) 
MAXAX2s 100 RATIO < 1 
RATIOS 2 DEGP 23 
DEGF 2 4 NAVEMB = 1 or 2 


DEPTHO/DEPTH= 0,5 67 DEP ITHQ=Z5 
NAVEMB = I or 2 


Grade = C (83) Grade = Cl (40) 
MAXAX2 Ss 100 RATIO? 4 
RATIO s 2 RATIO s 2 
DEGF 2 4 DEGF25 
DEPRAHQ/DEPIHS 0.5 NAVEMB 2 1 
DEPTHQ'> 25 


NAVEMB = 1 


Grade = D (197) Grade = DI (249) 
All run free which are All fixed depth which are 
not A, By. Or ©. not Al, BI, or Cl. 


The transition from Grade A to D is fairly clear. That is: 
A: Reliable hypocenter and My, 


B: Reliable hypocenter but perhaps less reliable m,, 


C: Reliable epicenter but perhaps less reliable depth as well as my, 
D: All others. 
The transition is more arbitrary and less clear for Al to DI. First note that DEGF = 5 
for AI, BI, CI, whereas it is DEGF 2 4 for A, B, C. In effect, this means we must 
have at least eight observations used for location in all cases. For AI, BI, we require 
RATIO < 1 rather than RATIO < 2 in the hope that tighter control on that parameter 
will compensate for not having DEPTHQ or MAXAX2. In case CI we include events 
with 1< RATIO < 2 in place of relaxing the depth control parameters which are not 
available. The transition in terms of NAVEMB is unchanged for I events. We should 
note that for many purposes events in grades Al, BI, and CI and even many in DI prob- 
ably have good depths. That is, shallow depths are probably good within a few tens of 
kilometers and deep events do get listed as deep. Unfortunately this is a subjective 
and unsubstantiated assertion. 

The definition of grades is somewhat arbitrary and the grades themselves are 
intended simply as aids to users of the event list. One can obtain some indication of 
the possible reliability of event parameters from the grades and for some research 
purposes one may wish to consider only events with high reliability to minimize the 
effects of errors in the event parameters. Unfortunately, for actual discrimination 
purposes, it would not be possible to consider only events with high reliability as de- 
fined by the grades. For example there are 30 D or DI events with listed m,, greater 


than 4.5, and 5 with m,, greater than 5.0. Even restricting interest to the Asian 


geographic region outlined on Figure 2 we find that of 40 events listed with m,, > 4.5 
there are three with grades D or DI. 

The 32 station list has also been graded to indicate the possible reliability of 
the 818 events which could be located with only that set of stations. Table [V summarizes 
the 32 station grades in comparison to the full ISM grades. In terms of numbers of 


events we see that the two lists are about the same in the low reliability categorics 


D and DI and the 32 station list has fewer events in the higher reliability categorics. 


TABLE IV 


Grade Distribution of ISM and 32 Station Lists 


Number of Number of 
Grade ISM Events 32 Station F vents 
A 306 Zar 
Al 71 77 
B 34 L2 
BI 16 7 
C 83 26 
‘@) 40 24 
D Oy 195 
i 249 240 


Fig. 2 Outline of geographic area defined by seismic regions 19, 
26, 2%s Ze, 29; 30, 41, 42, 48, 49. 


m, and Detection Statistics 

Figure 3 shows a cumulative histogram of the m,, values (AVI-MB) for the 859 
ISM events for which we were able to determine magnitudes. It also shows a histogram 
for events within the Asian region defined on Figure 2. In the selected Asian region 
there were 215 events with Mm), assigned and 33 with no mp, value determined. 

It is often assumed that the log of the number of cvents in any region and time 

interval will vary linearly with my One can then fit straight lines to histograms such 
as those on Figure 3 and estimate at what magnitude level the detection system secms to 


(9) 


start missing events. Basham and Anglin used preliminary ISM data for a region 
similar to that of Figure 2 and concluded that the cumulative 90°. detection threshold 
for that region was about 4.0. Since incremental thresholds are normally 0.2 to 0.3 


(6) 


units higher’ * they estimated the 90%, incremental threshold to be 4.3 for that region. 
The data of Figure 3 do not contradict the conclusion although this fitting of straight 
lines must be accepted with some scepticism as a method of estimating detection 
thresholds since it does use a theoretically unsubstantiated assumption about the linear 
variation of seismic activity with mM), 

Fitting the complete ISM cumulative occurrence of events to obtain a detection 
and location threshold is doubly bothersome. Not only must the straight-line assumption 
be accepted but it is even more difficult to find a straight region to fit with a straight 
line. However, if such a method is used, it is our estimate that the 90), cumulative 
threshold obtained is between 4.5 and 4.6 with the corresponding incremental threshold 
0,2 to 0.3 m,, units higher. 

Detection and location thresholds are discussed further on but first the ISM m 


hy 
(AVEMB) must be related to the PDE, LASA, and NORSAR bulletin valucs. Figure 4 
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Fig. 3. Cumulative histogram of complete ISM mM 
for events limited to Asian seismic regions 19, 26, 
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Fig. 4. 


shows histograms of the difference between the ISM mM, and the ms from these other 
sources. These data hold no special surprises. The scatter between PDE and ISM 

is less than in the other two cases since both PDE and ISM often average several 

station values to get their mp): The average deviation between ISM and PDE, LASA, 
NORSAR are 0.04, -0.05, and 0.14 units, respectively. Only the last, indicating NORSAR 
my, values are high by 0.1 to 0.2 units, might be considered significantly different 

from zero to be of any concern. In general there is scatter between these various 

mM, values but no large biases. 

Now that the relationship between these My values has been statistically charac- 
terized we can consider more data related to detection thresholds. Figure 5 shows a 
histogram of 1SM events which did not appear on the USGS PDE list. From this and 
Figure 3 we estimate that the PDE 90% cumulative threshold is about 5.0 with the 
incremental 0.2 to 0.3 units higher. This estimate assumes that the ISM missed no 
events above m, = 39.0 and that the scatter of my values can be ignored. It is con- 
sistent with what one would get using the fitting of straight-lines approach to PDE data 
itself. 

Figure 6 contains histograms of possible events which appeared in the LASA 
and NORSAR bulletins for the ISM time period but which could not be verified sufficiently 
by other station reports to be included in the actual ISM list. The counts for events 
placed in the region of Asia defined by Figure 2 are indicated by the horizontal lines. 
Of course these events are unverified, poorly located, and all mM), values are fora 
single station. Nevertheless if Figure 6 is accepted on face value it is in rough agree- 
ment with the detection thresholds discussed above for the ISM. First consider the 


Asia region only. A 90% incremental threshold of n= 4.2 or 4.3 is consistent with 


17 


70 


| 
50 


a qe 


30 


a iq 
40 | -—| ! ! 


Try 
OCULEELEALE EEE bp 7 


AVEMB 


NUMBER OF EVENTS 


Fig. 5. Histogram of m, (AVEMB) for ISM events not reported 
by USGS PDE, 


18 


NUMBER OF UNVERIFIED 


NUMBER OF UNVERIFIED 


. [18-2 -11684 | 


oe) 

lu 

= 1 4 
iz TOTAL = 35 

" i 

lu » 

s GP INDICATES NUMBER 
= IN REGION SHOWN 
_ 8 ON Fig.2 

=) 

>» 

2 | 

x 

= § 

Se) 

8 

O 

=< 

ee) 

WW TOTAL = 96 

r 

= 

< INDICATES NUMBER 
x IN REGION SHOWN 
= ON Fig.2 

uJ 

=4 

ad 

= 

ie) 

<q 

ep) 

= 6 

J 


_ 


Fig. 6. Histograms of unverified large-array bulletin events. 


(a) LASA my (LMB) for unverified LASA bulletin entries. (b) 
NORSAR my (NMB) for unverified NORSAR bulletin entries. 
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the fact that none of the potential events from that region on Figure 6 have assigned 
my, values greater than 4.2. If all the events on Figure 6 are real then these data are 
also not grossly inconsistent with the previously estimated worldwide 90%, threshold 
somewhere between 4.7 and 4.9. 

In general we believe that our data and methods do not justify assertions about 
detection thresholds which are any more precise than those we have made above. 

Figure 7 is one last mM, histogram. It shows the my, distribution of all ISM 
events with epicenters less than 95° and more than 30° from LASA. The shaded sec- 
tions represent the counts of the events which were not reported on the LASA daily 
bulletin. It can be seen that LASA reported a large fraction of these events and the un- 
reported ones, as a fraction of total events at any magnitude level, are distributed 
quite uniformly in magnitude. This would be the case if LASA were simply turned off 
for a fraction of the total ISM. In fact LASA, for various rcasons, was reported not 
operating for a time equivalent to 8.3%, of the ISM and the missing events represent 13.4) 
of the ISM events. An event by event check of all missing events showed that half of them 
were in the reported down times and the other were not. At the present time we have 
no satisfactory explanation for the missing events since, in our opinion, LASA should 


have detected most of them. 
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Fig. 7. Histogram of mp, (AVEMB) for ISM events located 30°-95" 
from LASA. Events not on the LASA bulletin are indicated by dark 
shading. 
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Location Comparisons and Evaluation 

Figure 8 is a comparison of the ISM epicenters of events with the locations 
given by other sources. These other sources are USGS PDE, LASA, NORSAR, and 
the 32-station list. Histograms of the differences in epicenters in degrees are shown. 
Agreement with PDE is excellent as expected. Most LASA and NORSAR epicenters 
are within 5° of the ISM location with nontrivial numbers of events out to 20° and 
occasional very large differences. The very large differences usually correspond to 
an incorrect phase assignment by the array in preparing the bulletin. The spread for 
the 32 station locations is more than for PDE but much less than for LASA or NORSAR. 
In fact the 32-station spread is canparable to that for the PDF when only PDF events 
are considered and becomes greater only when some of the more marginal events are 
added to the population. 

A comparison of ISM depths and those of the USGS PDE, is shown as a scatter 
diagram in Figure 9. That Figure contains points for every PDE event which was not 
assigned depth 33 by PDE and which was located by us without restraining depth. It 
is difficult to make conclusions from these data since the relationship of neither the 
PDE nor ISM depth to the true one is known. 


(3) 


As explained elsewhere’ “ our location programs find the origin time and 
hypocenter which minimize the goodness of fit parameter CHI2. This parameter is 
a Weighted sum of the squares of observation residuals. The weights are obtained by 
inverse scaling of residuals with assigned standard deviations before squaring and 
summing. An error ellipsoid can be defined which corresponds to a unit increase in 


CHI2 from the minimum value obtained in locating an event. Our hypocenter quality 


variables DEPTHQ and MAXAX2 are defined in terms of this ellipsoid as indicated in 


Ze 


NUMBER OF EVENTS 


t 
T 


- fi 
J r 
4 | 
i 
; } 
I 
( 

100+. | 
ml : 
zi n 
a H 
1 i 
4 : 
i | 
| — 
tl : 
{| ! 

10 | 


Boe lec 
ee Se ee 


= : aa 7 
i) 
bal 


Oo 


20 40 60 
PIDIST (deg) 


@ 
o 
oO 
Oo 


Fig. 8. Histograms of the distance (degrees) between the ISM 
epicenter and (a) USGS PDE epicenter, (b) LASA bulletin epi- 
center, (c) NORSAR bulletin epicenter, (d) 32 station ISM epi- 
center, for events with locations assigned by both sources 
involved. 
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Fig. 9. ISM depth (DEPTH) vs PDE depth (PDEPTH) for A, B, 
C, D grade ISM events with a PDE depth assigned other than 33 
km. All depths in kilometers. 
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Figure 10. It should be noted that DEPTHQ is not just the length of the line from the 
center of the ellipsoid to the surface of the ellipsoid in the vertical direction. For 
very long narrow error ellipsoids that might be quite misleading whereas DE} PTHQ 
would better reflect the true state of affairs. Similar comments apply to MAXAX2 
and in fact to the calculation of the hypocenter ellipsoid itself from the complete four- 
dimensional error ellipsoid which has a time as well as three space dimensions. 

Some experiments have been completed which give a more direct interpretation 
of the quality variables DEPTHQ and MAXAX2. One of these deals with DE PTHQ and 
the other with both DEPTHQ and MAXAX2. Results are described below. 

For the first experiment all A grade events which included depth phases in 
obtaining the ISM location were selected and relocated without using the depth phases. 
Of these 141 converged properly without any need for restraining depth tor the relocation. 
Since the previous locations with depth phases tended to pin the depth closely to that 
indicated by the depth phases we took the difference between the two depths as the 
actual error in depth for the relocated events. The absolute value of this error was 
weighted inversely by the depth quality of the relocated event and the histogram of 
Figure 11 generated. Accepting the assumption that the original depth was correct we 
see that the indicated depth of the relocated events was almost never in error by more 
than 2*DEPTHQ and that it was less than DEPTHQ for 75% of the relocated events. 

The second experiment also started with a set of A grade ISM events. for each 
event the basic data used to obtain locations were split into two approximatcly equal 
disjoint sets and two new locations calculated. The absolute value of differences in 
the two depths were calculated and normalized by the square root of the sum of the 


squares of the two DEPTHQs. Similarly the length of the vector between the two epi- 
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Fig. 10. Cross section through hypocenter error ellipsoid show- 
ing definitions of quality parameters DEPTHQ and MAXAX2. Sec- 
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Fig. 11. Histogram of the change in depth resulting from relocat- 
ing A grade events after discarding all depth phases. Each change 
is normalized by the relocation depth quality parameter. Total 
number of events used is 141. 
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centers was calculated and normalized by the square root of the sum of the squares of 
the two MAXAX2s for each event. The particular normalizations used assume that 
DEPTHQ and MAXAX2 are quantities conceptually similar to standard deviations and that 
they will combine in the same way as standard deviations when differences of the two 
independent estimates are taken. Although this is not strictly correct it was felt to be 

a satisfactory assumption for our preliminary calibration of Dk PTHQ and MAXAX2. 
Histograms of the normalized differences are shown on Figures 12 and 13. It can be 
seen that the normalized differences have about the same distribution as those of Figure 
I 

We would like to be able to claim that our calibrations of DF PTHQ and MAXAX2 
are absolute. Unfortunately we cannot really do this. Consider the last experiment for 
example and suppose an event is in an arc area with significant lateral inhomogeneities. 
Both of our relocations may have similar systematic errors and our calibration may 
thus be missleading. This question of absolute calibration is still essentially unresolved 
and will be the topic of future research. 

One final point to be made about ISM locations is that the reliability of assigned 
epicenters clearly deteriorates with decreasing body wave magnitude. Figure 14 shows 
a plot of MAXAX2 vs AVEMB for all ISM events for which at least three stations were 
used in calculating AVEMB and for which MAXAX2 could be calculated. Requiring threc 
Or more stations increases the reliability of AVEMB although errors in depth or location 
may still introduce nontrivial errors. It can be seen that the apparent reliability in- 
dicated by MAXAX2 begins to get considerably worst for significant numbers of events 
when AVEMB goes below about 4.7. This is not surprising if we recall that previous 
discussions indicated that this probably is about the 90°, worldwide incremental thres- 


hold for the ISM. 
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Fig. 12. Histogram of the normalized absolute value of the differ- 
ence between depths obtained by locating the same event using dis- 
joint sets of observations. Normalizing factor for each is the 
square root of the sum of the squares of the two DEPTHQ quality 
factors obtained for the event. One normalized depth error was 
off scabe at 5. 3 
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Fig. 13. Histogram of the normalized difference in epicenters 
obtained by locating the same events with disjoint sets of observa- 
tions. Normalizing factor for each is the square root of the sum 
of the squares of the two MAXAX2 quality variables obtained for 
each event. 
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Fig. 14. Epicenter quality parameter (MAXAX2) vs mp (AV EMB) 
for all ISM events with m, determined from an average of three 
Or more stations. 
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©; DATE 
O, TIME: 
LAT: 
LONG: 
DEPTH: 
AVE MB: 
NAVEMB: 
MAXMB: 
MIN MB: 


NEXMB: 


DEGF: 


CRir2: 


RATIO: 


MAXAX2: 


DEPTHQ: 


GRADE: 


REGNUM: 


SE Lous Ga 


NAT: 


APPENDIX I 

Parameter Definitions 
Date (Month/Day/Year) 
ISM origin time (Hours:Minutes:Seconds) 
ISM Latitude (Degrees with + for North) 
ISM Longitude (Degrees with + for East) 
ISM Depth (Kilometers) 
ISM Mm), (Average of single station mM, values) 
Number of stations averaged to get AVEMB 
Maximum mM), value calculated for any station 
Minimum My, value calculated for any station 
Number of station my, values excluded in calculating AVEMB. Station 
excluded if its my value is more than 0.6 units from AVEMB 
Degrees of freedom for the ISM location and origin time estimate. 
Sum of the square of residuals of observations from ISM origin time 
and location. Residuals weighted inversely with assigned variance. 
CHI2/DEGF 
Length of the major semiaxis of the ISM epicenter error ellipse (kilo- 
meters). See Figure 10 
Half length of depth error bar (kilometers). See Figure 10. 
ISM event grade assigned as defined in Table 
Geographical region number of ISM event 
Seismnic region number of ISM event 


Number of station arrival times associated with ISM event. 


NV: Number of array velocity observations associated with ISM event 


NATLOC: Number of station arrival times used in locating ISM event 

NV LOC: Number of array velocity observations used in locating ISM events. 

POTIME, PLAT, PLONG, PDEPTH, PMB: USGS PDE values corresponding to O. TIME, 
LAT, LONG, DEPTH, AVEMB. 

PIDIST : Distance from ISM epicenter to USGS PDE epicenter (Degrees) 

LOTIME, LLAT, LONG, LDEPTH, LMB: LASA bulletin values corresponding to O. TIME, 
LAT, LONG, DEPTH; AY.EMB, 

LIDIST: Distance from ISM epicenter to LASA bulletin epicenter (Degrees) 

NOTIME,NLAT,NLONG,NDEPTH, NMB: NORSAR bulletin values corresponding to 
O. TIME, LAT, LONG, DEPTH, AVEMB 

NIDIST: Distance from ISM epicenter to NORSAR bulletin epicenter (Degrees) 

XOTIME. . .XNVLOC: 32-station values for similarly named ISM variables without 


the leading "X”"’. 


ISM EVENT LIST 
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Month —1SM) has been compiled using short-period data from a large number of seismic stations and 
arrays. Hypocenters and other parameters have been obtained for every event and are presented 
along with comparisons with other event lists for the same time period, Also included are the re- 
sults of experiments designed to calibrate certain parameters which have been selected to indicate 
the accuracy of estimated depths and epicenters. The ISM list contains 996 verified events. 
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